Fermi Edge Singularities: Bound-states and Finite Size Effects 
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Fermi edge adsorption singularities (FES) are studied using a combination of conformal field the- 
ory (CFT), an exact sum rule and numerical work on a tight binding model which is shown to exhibit 
remarkable simplifying features. The relationship between FES and the Anderson orthogonality ex- 
ponent is established in great generality, using CFT, including the case where the core hole potential 
produces a boundstate. Universal results on the adsorption intensity in a finite sized sample are 
obtained. Various predictions are checked numerically and the evolution of the adsorption intensity 
with electron density is studied. 

I. INTRODUCTION 



5_H 

O |' A theoretical understanding of the Fermi edge singularity (FES) in X-ray a ds<WjBt»a* 

^ ' als and of the related Anderson orthogonality catastrophe dates back to the 1960's.lifEiU Q t 

Nonetheless, this remains an active area of research today, in part because of theoretical and ex- 
perimental worte^H-i |y(j j tnm|i q f-Hre4uced dimensionality where strong correlation effects may play an 
important roleJiirll3'tjllliEj't3E3'ta In particular pinew and very general theoretical approach has 
been developed, based on conformal field theory.t3 

Recent experiments have studie*Uie*>tical absorption threshold singularities in modulated, doped 



semi-conductor layered structures. LjliS These behave more or less like two dimensional metals with a 



1 continuously tunable electron density. The potential produced by the valence band hole is expected 

' to produce a boundstate (exciton) due to the well-known theorem that an attractive potential in 

, two dimensions always has a boundstate. Many interesting issues are raised by these experiments 

which have not yet been adequately addressed theoretically. In particular, it is possible to study 
t**"* \ the behaviour of the threshold singularities as the electron density goes to 0. In some cases an 

ON ■ additional threshold corresponding to a negatively charged exciton (two conduction electrons bound 

to a valence hole) is observed. Strong correlation effects may play a crucial role in this case and the 

S usual Fermi liquid approach may need to be modified. 
The present work addresses several issues in this field, within the usual Fermi liquid framework. 
The case of a core potential producing a boundstate is considered, as is the behaviour of the threshold 
singularities as a function of electron density. Furthermore, we consider the nature of the adsorp- 
tion intensity for a simple model of a finite sample. The techniques employed in this paper are a 



a 

^ | combination of conformal field theory methods, an exact sum rule and numerical work on a tight 

binding model which exhibits remarkable simplifying features making it feasible to study very large 
systems. 

, The conformal field theory method of Ref. ( [12|) is extended to the case where there is a bound- 

state. A very simple and general proof of the exact correspondence between the FES exponent and 
the orthogonality exponent is given. The adsorption intensity near threshold, for a finite system is 
shown to have a simple universal form, using conformal field theory. An exact sum rule is introduced 
which determines the ratio of adsorption intensities for the two cases where the core bound state is 
empty or filled. 

Numerical work on the simple tight binding model is used to check the validity of Anderson's 
formula relating the orthogonality exponent to the phase shift at the Fermi surface, the Nozieres-de 
Dominicis, Combescot-Nozieres formulasB'H for FES exponents including the boundstate case, and 
the formulas newly derived here for the adsorption intensity in a finite system. In addition this 
numerical work, together with the sum rule, is used to study the behaviour of the adoption intensity 
as a function of density. 

Some of the new conformal field theory results were briefly described in Ref. ( |l9| ). 



II. CONFORMAL FIELD THEORY APPROACH 



An approach based on boundary conformal field the^ytjt-j provides a unified view of the problem. 
As usually, we start from the simplest possible modeln, 
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Here the band of spinless, non-interacting electrons (c',c) is scattered by the core hole potential 
Vk,k'- The core hole (b\b) is dispersion-less. It is created and annihilated instantaneously (by 
irradiation), and the reaction of the band electrons on this instantaneous perturbation constitutes 
the FES and Anderson orthogonality effects. 

In the low-energy limit, the system is mapped onto the (l + l)-dimensional Dirac Fermions defined 
on a ray r > with a scattering potential V at the origin. (This can be done by assuming spherical 
symmetry of e k and V ki y and considering only s-wave scattering; generalizations to other cases are 
straightforward.) The following discussion applies, with minor modifications, to either an s-wave 
projected 3D (3 dimensional) problem or to a problem defined a priori in ID. We will henceforth 
generally consider the ID case. We could consider, for example, a ID tight-binding model, defined 
on the positive half-line, with free boundary conditions and a potential localized near x = 0. See 
Appendix A for a detailed discussion of this model. The corresponding boundary condition in the 
low energy Dirac theory is: 

Vi(0) = Vfl(0). 

The role of the scattering potential in the theory here is to impose an effective boundary condition 
on the low energy degrees of freedom, relating the left and right movers: 

Vfl(0) = e 2a(fcF) <M0), 

Here 8(kp) is the phase shift at the Fermi surface; fc-dependence of the actual phase shift is irrelevant 
at low energies. The action of the hole creation operator, b\ thus reduces to that of a primary 
boundary condition changing operator, 0. The Green's function (hole propagator) of this operator 
in a half-plane, z = r + it, r > 0, is 

G(n -r 2 ) ee {b{T,)b\T 2 )) = (A; 0| 0(n)0(75) |A;0> = ^ } ^ . 

Here x is the scaling dimension of 0, and |A;0) is the ground state of the infinite system (filled 
Fermi sea) without scattering potential. Physically, the Green's function is directly related to the 
absorption intensity in the case of photoemission, 

J(w) oc J dte l{uJ -" a)t {b(t)b\Q)) oc (w - w )"°, 

where ljq is the threshold frequency. Evidently, the (FES-) exponent a and the scaling dimension 
are related via 

7 ee 1 — a — 2x. 

Looking for finite size effects, we conformally map the half-plane onto the strip, I > r > Using 
the transformation z = le ww ^ 1 . In a bosonic system this automatically gives the same boundary 
condition at and I. However, for fermions it gives: 

<M0) = tM0), i>L{i) = -^ R {i). 

This follows because the fermion fields transform as: 

TpL — » (dz/dw) 1 ^ 2 '4'L 
tp R — » (dz* /dw*) 1/2 i) R . 

At w — x + il, 

(dz/dw) 1/2 /(dz*/dw*) 1/2 = -1. 

This transformed problem corresponds to considering a ID model defined on a finite line, < 
x < I, with the impurity potential near x = and an appropriate boundary condition at x = I. 
Alternatively, the 3D s-wave projected system is now defined inside a finite sphere of radius I with 
an appropriate boundary condition on the surface of the sphere. For a discussion of this boundary 
condition and more details, see Appendix A. 
We find for the Green's function on the strip 



(AA;0\O(u 1 )O(u 2 )\AA;0) = — (2.10) 

(# sinh i («i-«a)) 

A^O) being the unpertur bed ground state of the system of length i, with the "same" boundary 



condition, A, g iven by Eq. (2.3) at both ends. 

In Eq. (2.1C) we can either Taylor expand sinh in the limit 7r(ui — U2)jl <C 1, or insert a complete 
set of states \AB;m) (eigenstates of the system with the scattering potential - boundary condition 
"B" - present), and obtain the relation 



(t 



e ' l + 2:re ' H i '-e * 



2:e(2:e + \)(2x + 2) ^^1-^2) 
H e ; 

D 
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= ^ \( AA > °l \AB;m)\ 2 e -l< B -< A ll^-^\ (2.11) 



If for the operator O the first nonvanishing matrix element is with the ground state of the 
perturbed system, |AB;0), then for the overlap of the two ground states (Anderson orthogonality 
catastrophe) 

|<AA;0|C?|^i3;0>l = (j)"- (2-12) 

The Anderson orthogonality exponent coincides with the scaling dimension x; on the other h; 
is given by the 0(1/?)- contribution to the ground state energy shift due to the perturbation,!! 

x = -[E^ B -E$ A \. (2.13) 

Here it is being implicitly assumed that this energy difference consists of a term of 0(1/1) only, as 
would follow from conformal invariance. As discussed in Appendix A there will in general also be a 
term of O(l) which must be subtracted. 

A simple way of determining the Fermi edge exponent 7 and the orthogonality exponent x, is 
thus to calculate the \jl finite-size correction to the difference in groundstate energies of the system 
with and without the scattering potential. The term of O(l) is non- universal, (cut-off dependent) 
while the higher order terms contain the corrections from various irrelevant operators. On the other 
hand, we expect the term of 0(1/1), which is determined only by the immediate vicinity of the 
Fermi surface, to be universal and to give the desired FES and orthogonality psponents. In fact, 
this result remains true including bulk Coulomb interactions in one dimensions Thus calculation 
of the 0(1/1) term in the groundstate energy difference gives a very simple way of determining the 
FES and orthogonality exponents in great generality. This calculation is spelled out in detail, for a 
one-dimensional tight-binding model, in App. A. The conclusion is: 
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While various derivations of this result, both for the FES exponent and for the orthogonality 
exponent have been given before, this one has certain distinct advantages. The original derivation 
of the orthogonality exponent by Anderson made a variety of approximations, including Taylor 
expanding certain quantities in powers of the phase shift. The derivation of the FES exponent 
in[ |^| also initially assumed a small 5(/cf) and then argued for the generality of the result by 
some fairly subtle consistency arguments. The bosonization derivations start with a bosonized 
Hamiltonian written in terms of <5(/cf) whereas a naive bosonization in fact only picks up the Born 
approximation to S, linear in the scattering potential. It is expected that eliminating the high energy 
modes somehow renormalizes this parameter in the bosonized Hamiltonian, turning it into the true 
phase shift. Once the assumption of conformal invariance is made, it is very straightforward to 
demonstrate that it is precisely the phase shift at the Fermi surface which enters the exponents, 
by an explicit calculation of the groundstate energy, as given in App. A. We note that once the 
bosonized Hamiltonian is assumed, the results for the strip can be obtained by a mode expansion 
of the boson field. This, of course, gives the same result obtained more simply by the conformal 
transformation. 

Another advantage of this somewhat abstract approach to the problem is that it can be immedi- 
ately generalized to the case where the core potential creates a boundstate. The Green's function 
can then be presented as a sum of two terms: 



(2.14) 
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G(u) = G e (u) + G f (u) 



\{AA;0\Q\AB;m;e)\ 2 e' 



+ y2\(AA;0\O\AB;n; f)\ 2 e 



where the first sum is taken over all states \AB;m;e) where the boundstate is empty, and the 
second over the states \AB;n;f) where it is filled. These two terms give rise to two peaks in the 
absorption rate, separated by the binding energy, Aco — \Eb\- Introducing the occupation number 
of the boundstate, fiB, 

G{u) =< [&(u)n s (u)][ns(0)6 t (0)] > + < [&(it)(l - nu(u))][(l - n s (0))6 t (0)] > . 

In the long time limit we may calculate each of these terms separately using the boundary con- 
formal field theory approach. The boundstate is associated with a finite binding energy, Eb and 
an exponentially decaying wave-function. Thus it has no direct effect on the 0(1/1) terms in the 
energies. Therefore we expect the above formulas to apply immediately for the first threshold where 
the boundstate is filled, with the exponent: 



5(e F ) 



When the boundstate is empty, the only change in the low energy physics is that one additional 
electron is raised to the first unoccupied state above the Fermi surface. This has wave-vector: 



k — k F + j 



We may regard b[l — ns] as a different boundary condition changing operator which creates one 
additional low energy electron, in addition to producing the new boundary condition of Eq. ( ^.3| ) . 
The 0(1/1) term in the "groundstate" energy difference in the case of the empty boundstate is, from 
App. A: 



v> v - „ n 1 J ( eF ) 

tb — Hjq — VF — - 

Thus the orthogonality exponent giving the overlap between the unperturbed groundstate and the 
"groundstate" with the boundstate empty is: 



5(e F ) 



The FES exponents for the two thresholds a/ and a e are given by: 

Q.f = 1 — 2xf, a e = 1 — 2x R . 

By merely computing the groundstate energy, rather than attempting to compute the exponents 
directly, we have finessed the problem of attempting to boscmize the theoryr-with the boundstate. We 
note that these results agree with Combescot and Nozieresfj and Hopfieldu. The present derivation 
seems quite closely related to the observation of Hopfield that the FES exponent measures the 
amount of charge pulled in from infinity by the boundstate. 

The conformal mapping from the plane to the strip establishes in a simple way the relationship 
between the FES exponents of the infinite system and the orthogonality exponents and energies of 
the finite system. In fact this mapping provides considerably more information. Let's imagine a 
rather artificial situation where a core hole is instantaneously created at the end of a finite one- 
dimensional system. (Equivalently we could consider an artificial situation where it is created at 
the centre of a finite sphere.) In this case the adsorption intensity, I(uS) of Eq. (2.5) becomes 
a series of 5-function peaks, as we see from Fourier transforming Eq. (2.11). The peaks occur at 
Et B -E$ A = n(x+m)/l, the energies of excited states of the perturbed system measured from 
the unperturbed gro und s tate energy Eq A . [The neglected term of O(l) just shifts the threshold by 
wo-] It follows from (2.11) that the ratio of the mth peak to the zeroth one depends only on x: 



\(AA;0\O\AB;rn)\ 2 _ 2x{2x + l)(2x + 2) . . . (2x + m - 1) 
\{AA;0\O\AB;0}\ 2 ~ m\ 



the O^ 1 peak, considered by Anderson, given by Eq. ( 2.12 ), 
independent pure universal numbers, determined only by <5(£;f). 

In fact, each of these peaks corresponds, in general, to several different states with energies that are 
degenerate, to 0(1/1). These are simply multiple particle-hole excitations of the free fermion system, 
with a dispersion relation which is linearized and a phase shift which is assumed fc-independent, near 
the Fermi surface. A general particle-hole excitation may be constructed by first raising n m electrons 
m levels, then raising n m _i electrons m — 1 levels, etc. The energy relative to the groundstate of 
this state is: 



Note that, from Eq. (2.11), each of these peak intensities scales with length the same way as does 

The ratios of peak intensities are length 



7T \ 

= 7 2^, 



n p p. 



p=i 



The first excited state has m = 1, n p = for p > 1. The next degenerate pair of states have ni = 2 
or ri2 = 1 (with the other n p =0). The third set of excited states is threefold degenerate with 
ni = 3 or m — ri2 = 1 or n-j — 1 (and the other n v — in all cases). Corrections to the linear 
dispersion relation and variation of S(k) near Icf will split these energies by amounts of 0(l/l 2 ). 
The simple prediction obtained here from a conformal transformation does not give the amplitude 
of each peak separately, but only the sum of amplitudes of all peaks at a given energy, where energy 
differences of 0(l/l 2 ) are ignored. 

This new finite size result interpolates, in a sense, b etwee n the orthogonality exponent and the 
FES exponent. Considering the large m limit of Eq. ( 2.22 ) we find that the intensity decays as 



-(1-2*) 



recovering the FES exponent. 



This result applies immediately to the peaks corresponding to the boundstat e bei ng fill ed or 
empty, provided that the appropriate orthogonality exponents, x; and x e of Eq. (2.17) and (2.20) 
are used. 



(2.23) 



III. ONE-DIMENSIONAL TIGHT-BINDING MODEL 



An evident discrete counterpart to the system (2.1) in its one- dimensional version is the system 



of spinless fermions on a finite ID chain with nearest-neighbour hopping, free boundary conditions 
and an impurity potential which can be switched on/off at the first site: 

H = H + b' f bHi; 

1-2 

Ho = -t^2 (VIV^+i + ; 
1=1 

Hi = -Vil>llK. (3.1) 

Here we choose V to be positive in case of an attractive core potential, and also choose t > 0. 
This model is very amenable to large scale numerical work with a minimum of effort. Not only can 
the single particle energies and wave-functions for finite I be found exactly in a simple form, but, 
more remarkably, the overlaps of the single particle wave-functions corresponding to different values 
of the potential, V, obey an exact factorization. This enormously simplifies the calculation of the 
overlap of the many-particle, Bloch determinant, wavef unctions. This can then be expressed by the 
Cauchy determinant formpla, used as an approximation by Anderson in his classic paper on the 
orthogonality catastrophe£3 The many-body overlaps are easily evaluated numerically and in some 
limits analytically (e.g. in the narrow band limit, — — » 0). As a result, chains of length up to few 
thousand sites are easily handled on a workstation using the standard "Mathematica" package. The 
crucial factorizability of one-particle overlaps disappears under any other position of the scattering 
potential. 

The details are given in the Appendix B; here we just summarize a few salient features, beginning 
with the infinite I limit. This model has a band of eigenstates with exact wave-functions: 

9j ocsin[fcj+«5(fc)], j =1,2,3,... (3.2) 

The dispersion relation is: 

e(k) = -2tcosfc. (3.3) 
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The phase shift is given by: 



S — arctan 



sin k 
t/V — cos k 



Note that at the bottom of the band, k — > 0, 8 — *■ for V < t, when there is no boundstate, but 
S — > n, for V > t when there is a boundstate, as required by Levinson's theorem. As k ranges over 
the whole band, from to tt, < S(k) < n/2 for V < t, and < 5(k) < tt for V > t, as shown in 
Fig. (|l]). For V ^> t, S(k) ~ -k — k. There is one boundstate, if V > t only, with: 

E B = -(V + t 2 /V). 

Note that this approaches the bottom of the band, — 2t, as V — » where the boundstate disappears. 



(3.4) 



(3.5) 



CO 




0.4 0.6 

FIG. 1. Phase shift in the ID tight-binding model vs. density v = lim;.,^ N/(l - 1), at t/V = 0, 0.05.0.01, . . . , 2, 3, 4, 5. The 
dotted line t/V = 1 separates regions with and without boundstate. 



The exact boundstate wavefunction is: 

*f oce~ Kj , 

with 

K = \a(V/t). 
For finite I, there is a set of wavefunctions: 

oc sinfc„(j - /), 
with the allowed wave-vectors determined by: 

sin k(l — 1) t 
sin kl V 



For t/V > 1 — l/l there are Z — 1 solutions of Eq. (3.9), which we label 



fei, ki, . . . ki-i. 

For t/V < 1 — l/l there are only I — 2 such wavefuctions, which we label 2, 3, ... (I — 1) and an 
additional wavefunction: 

i&j oc sinh «(j — I), 



(3.6) 
(3.7) 
(3.8) 

(3.9) 
(3.10) 

(3.11) 



with k the solution of: 



G 



sinh k(1 — 1) t 
sinh td V 

In the case V = 0, the wave-vectors are: 

fc„ = 7rn/Z, n — 1, 2, 3, . . . (I — 1). 



From Eq. (3.9) we see that, in the limit V/t — > oo, the wave-vectors are: 

kn = (n- 1)tt/(Z - 1), n = 2, 3, 4, ...(«- 1), 

corresponding to a chain of (I — 2) sites and a free boundary condition. The lowest wavefunction 
(boundstate) becomes localized at j — 1 in this limit with eigenvalue —V. 

The overlaps between the unperturbed, ip, and perturbed, tp, one-particle states assume the 
special form, reminiscent of the first-order perturbation theory, but actually exact (Appendix B): 

< $ > = _yC{k m )C(k n ) 



e(k m ) - e(k n ) ' 



where 



C(fc) = 
C s = 



e B — e(fc„) 
v / 2sin(Z - l)k 



(l — V) — sin(Z — l)fc cos Ik/ sin fc 
v / 2sinh(/ - 1) 
•^/sinh(Z — 1)k cosh Z/t/ sinh k — (Z — 1) 



IV. CALCULATION OF HOLE PROPAGATOR. COMPARISON TO THE CFT PREDICTIONS 



The Green's function G(u) is determined by the set of matrix elements, 
\(AB;rn\O\AA;0}\ (4> m;e ,/| *o> = C e °,/ 

between (perturbed and unperturbed) many-body states of the system and corresponding excitation 
energies, Ae m;e j (we note explicitly whether the boundstate if present is occupied (/) or empty (e), 
and label by m' the appropriate excited states of the band). 

When there are TV spinless noninteracting electrons in the system, the many-particle wave function 
is a N x N Slater determinant 



$ = 



1 



det(*J 



Here ^ ( j nb ^ is an appropriate one-particle eigenfunction of the state rib, taken at the coordinate of 
the ath particle, a,b = 1,2, ... N. 

The overlap of two such states is a determinant 

($,$) = i_^...^det(^™ ) *)det(*<7 ) ) EEdet((* (m) ,* (n) )) . 



The pepiarkable form of the one-particle ov erlap s (3.15, 3~16| ) allows us to apply the Cauchy 
formulaE2l in order to calculate the determinant ( O ) : 



det 



V a m + b n J 



TT ^ (a m — a n )T\ . [bn 



fen) 



]lm,«( a ™+ 6 ») 



ieresQi, 



Unlike the situation considered in the papers by Anderson and Combescot and Nozierestj'U, m our 
model the special form of the overlaps is an exact result, and not the consequence of the linearization 
of the dispersion law close to Fermi surface. 

We begin with calculating Cf°, the overlap of the grou nd st ates of the system with and without 
the core potential, which yields the Anderson exponent (2.12). In the corresponding determinant 
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(4.3) both "old" and "new" indices run from n, m = 1 to n, m = N. (That is, in the "new" state 
there is one bound electron and (N — 1) electrons in the band, and no e — h pairs). The other 
interesting overlap, C°°, corresponds to the situation when the bound state in the "new" system is 
empty, and all N electrons are in the band, occupying the lowest lying states (still no e — h pairs): 
n = 1, . . . , N, but m = 2, . . . , N + 1. It should yield the "empty boundstate" Anderson exponent, 
which according to general considerations!! should be (l — ^) , as distinct from the "filled" value 



2 

) . Using the Cauchy formula, we find 



C 00 = { _y fGB -Q c( - m) -Q O^ IC^! ^)Un>nU^' *J (4B) 



m=2 n—1 

N+l N T1 N + 1 l~ - \T~f N 



C e ={-V) [[ C{ Km ) J_J_C(k„) +1 — , (4.6) 

m = 2 n=l llm=2 lln = lV £ ™ £ ™ ) 



and for their ratio: 

_ C"" _ C(KN+l) 61 — 61 TT gjV+1 — 6m TT Cm — £l ^ ^ 

J m — 9 ™ — 9 



m — 2 m=2 



Expression ( [f.Tj ) is easily calculated for pretty large systems, since it involves the number of opera- 
tions only of order I. The approximate values for the energies E m can be accurately calculated as a 
perturbation series in t/V <C 1 (see Appendix B). If the above predictions are valid, then the ratio 
should depend on the system size as 

Roc (I- 1) 2 U ) 2 I 1 =(l-i)—s. (4.8) 
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FIG. 2. (a)Logarithm of ratio \R\ = |C°° /C° f °\ as a function of (even) chain length (I - 1) for t/V = 0.1 and 
v = 1/2; 2/5; 1/4; 1/10 and 1/100; the lines are best fits to points with (I - 1) > 10 by In \R\ = A + (1n(l - 1); (b) Ex- 
ponent C(u). The curve is ^-H - 1 ) _ i ( c )The same as (a) for t/V = 0.01 and v = 1/2; 1/5; 1/11 and 1/101 (in the latter 
three cases the length (I — 1) was chosen to be odd); the lines are best fits to points with (I — 1) > 10 by In \R\ = A + £ ln(l — 1); 
(d) Exponent ((v). The curve is =°- 01 ) _ !_ 
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Here 5f is the phase shift at the Fermi surface, k F = kv. 

The results are shown in Fig.^ (t/V — 0.1 and 0.01) for (/ — 1) < 4000. As is clear from the figures, 
the ratio as a function of I at fixed density v — behaves as 

R(l;v)<x(l-l)« u \ 

where indeed (the least squares best fit parameters are shown in figure captions): 

5 F (u;± =0.01) 1 



CM 



for ~ I _ Vj as seen in Fig.l. 

Calculation of the coefficients Cf° e themselves is more time consuming. We calculated C®° and 
Cf in the limit t/V -> for (/ - 1) < 400. The results are shown in Fig|| The matrix elements 
dependence on size and density is accurately described by 

_i / i'um V 

c° e ° f (i,u)=A e , f (u)(i-i) n - ; , 

thus confirming the validity of the original Anderson's result0 in the case of binding core potential 
and arbitrary electronic density. 



(4.9) 
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(4.11) 
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(a)Logarithm of the coefficient \C^°\ 
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(empty boundstate) as a function of chain length for t/V —* and 
v = 1/1000; 1/500; 1/100; 1/20; 1/7; 1/5; 1/4; 1/3; 1/2; 3/5 and 4/5; the lines are best fits to points with (I - 1) > 10 by 
C°°(l,u) = A e (u)N~^" M ; (b) The same for |C°°| (filled boundstate). (c,d) Anderson exponents /3 e (c) and f3 f (d). The curves 

are (3 e {u) = and ^/(i/) = il^ll. (In the limit t/V = 0, v = 1 - 5/tt.) 



Turning to the contributions of the excited states, CS m , we must keep in mind that in the model, 



the degeneracy of energy levels in the CFT formula ( 2.1?: ) is lifted. Therefore when checking Eq 



iE' c °l\ 

(2.22), we calculate the ratios J — — f^— , where the sum is extended over the excited states which 



would be degenerate in the case of linear dispersion relation. These clusters of nearly degenerate 



states are clearly seen in Fig. (|4j), where we plot log f(E) = log y"]_ 



<E | 



vs. log E. In fact, 



9 



not all excited states are included in Fig. (|4|), but only those corresponding to a single particle- 
hole excitation. In the large I limit, all excited states with energies nVpm/l with m = 1,2 or 3, 
discussed in Sec. II, correspond to single particle-hole excitations. However, multi particle-hole 
excitations begin to appear at m = 4. We might try to extract the FES exponent from the slope 
at low energies and large I, assuming that the neglected multi particle-hole excitations don't make 
too large a contribution close to the threshold. Assuming I ~ La~ a in this region, then for small E 

f(E) ~ E 1 -" = E 1 . The tangent of the curve should yield the FES exponent-related 7 = , 

but actually it is significantly smaller. 



0.4 



0.3 



LU 

' S 0.2 



0.1 








In E 

FIG. 4. The normalized integral of intensity, f(E) = , 



I(uj)du>, vs. energy, for the case of empty 
168. The best fits yield 7 = 0.037; 0.047; 0.054 



■iE m <E | C 00| ; 

boundstate and t/V — ► 0. Density v — 6/21; system size (/ — 1) = 21; 42; 84; 
and 0.060 resp. The CFT value is 0.082 (uppermost line). The clusters of almost degenerate excited states are clearly seen (in 
one-pair approximation the n-th excited state is n-fold degenerate). The dotted line shows the direct diagonalization results by 
Eder and Sawatzky (Ref.( 0)), for t = 1, V = 32, N =6,(1-1) = 21; 7 = 0.021. Due to normalization to | C°° | 2 the curves are 
offset in vertical direction; as the system size grows, the plateau heights are reaching the univers al valu es predicted by CFT for 
the ratios of subsequent excited peak amplitudes to the lowest energy peak amplitude (see Eq.(2.22) and Fig.^). Inset: FES 
exponent vs. inverse chain length. The cubic extrapolation to infinite system size 7(0) = 0.068 still falls short of the CFT 
value. 



The calculations of the intensities of the first few peaks ( t/V = 0.1; Fig.g) shows that finite-size 
corrections to these amplitudes are significant for I ~ 400. Thus we shouldn't expect to be able to 
obtain the FES exponent reliably this way even if the neglect of multi particle-hole excitations was 
valid. 

In Fig. (|H|) and the Table we show the ratio of the total amplitude of all peaks at excitation 
energy VFirm/l, to the amplitude of the lowest peak (m = 0) for m = 1,2,3 for both filled and 
empty boundstate, JJ^, following the discussion in Sec. II. t/V = .1 and the ma ximu m length 
considered was / = 400. The predictions for R%f from conformal field theory in Eq. (2.22) are also 
shown in the figures and table. The finite-size corrections to these ratios are quite large, but upon 
extrapolating in 1 /I we obtain good agreement with the CFT predictions in all six cases. 



10 



CD 




CO 




"go 


"O 




CD 




O 


0.1 


X 


CD 








to 








CD 
-C 


0.8 






O 


0.7 








0.6 


CD 




s 


0.5 


CD 


> 
CO 


0.4 


CD 
CC 


0.3 




0.2 




0.1 





(a) 


^0.2 












0.5 




0.6 , ' 





CD 
3 

GO CD 



0.25 
0.2 



a 
o 

_Q 

>, 
Q_ 

E 

LU 



2 0.15 

T3 

£ 0.1 

'o 

« 0.05 



\v=.0.1 
\ 


(b) 


\ . 
^. .0.2 \ , 




^0.3 ^* " . 




0.4 ; * 




■0.5 ' 









.02 



0.02 


0.04 0.06 
1/(1-1) 


0.08 


0.1 


.. v=.0.6 






(c) 


S 0.5 y 








-. 0.4 ' • . 








. 0.3 . 








.02 








0.1 

0.02 


0.04 0.06 


0.08 


0.1 



0.02 



0.04 0.06 
1/(1-1) 



0.08 



CD 

to 

00 
T3 
c 

=5 

o 
E 

LU 




CD 

to 
T3 
C 
Z5 

o 

Q. 

E 

LU 




0.1 



0.04 0.06 
1/(1-1) 

FIG. 5. (a) Relative weight of the first excited state peak in the hole propagator for different densities 7V/(Z — 1) = 0.1(0.1)0.6; 
t/V = 0.1. The horizontal lines mark the CFT predictions (Eq.( |2,2^ )), curves are the least squares' fits of the third-order in 
1/(1 — 1), calculated from the four points with smallest 1/(1 — 1). Notice that these curves are in a good agreement with the 
data in a larger range of 1/(1 — 1). (b),(c) The same for the second and third excited states. 



The finite size corrections, at a fixed value of I, are roughly proportional to 




, where the 



effective phasspphift in the presence of bound state is Sf — 8f if the boundstate is filled, and tt — Sf, 
if it is emptyau. Thus CFT predictions in a finite system with filled boundstate are more accurate 
for low density than for high density, and the opposite for empty boundstate. 
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V. SUM RULES AND OPEN QUESTIONS 

The Fourier transform of the core hole Green's function, proportional to the photoemission in- 
tensity, can be written: 

I{u, v) = ! (0|n) \ 2 5(E n -E - u>), 

n 

where |n), E n label all states of the system with the core hole potential turned on and |0), Eq refers 
to the groundstate without the core hole potential. Here we make explicit the fact that I depends on 
the electron density (i.e. fci?). The states |n) can be classified according to whether the boundstate 
is filled or empty and accordingly we may decompose 1(d): 

I(u>,u) = If(u,i/) +ie(w, v). 

Near a threshold, I(lo, v) takes the form: 

I(u,u) =F{v){io-uj )- a{v) . 

So far, we have focussed on the value of the FES exponent, a(v). In this section we would also like to 
consider the dependence of the amplitude factor, F on the density and also the behaviour of 7(a), v) 
away from the threshold. While most of this behaviour is clearly non-universal, one might expect 
certain universal features to immerge in tbe-ilimit v — > 0. Our interest in this limit is motivated by 
the experiments on doped semiconductors.!^ 

To begin with we point out the existence of two sum rules. These apply very generally to FES 
problems in arbitrary dimensions, without any particular assumptions about spherical symmetry or 
about the location of the core potential in the one-dimensional case. It follows from completeness 
of the states |n) that I(u)) obeys the sum rule: 




du>I(uj, v) = 1. 



This implies that the integrated intensities from the states with filled or empty boundstate obey: 

// + Ie = 1. 

Another useful sum rule can be derived by writing 7/(cj, v) in terms of the projection operator, hs 
onto states in which the boundstate is occupied: 

J dull f (u), v) — ' I (0|ns |n) | 2 = (0|ns|0) . 

n 

We may write 

where i/)b annihilates the boundstate electron. ipB can be expressed in terms of the operator ijjj 
which annihilates an electron at site j and the boundstate wave-function, 'iff : 

3 

Thus: 

/ / = ^*f*f <o|v4^|o). 

In a d-dimensional continuum formulation this becomes: 

I } = J d d rd d r^ B {^ B {r)(Q\^\r)^{r)\Q). 

At this point, it is convenient to Fourier expand the position-space annihilation operators. The form 
of this expansion depends somewhat on the boundary conditions (in the system without the core 
potential) . In dimension d > 1 we generally consider a translationally system in which case: 
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** = j Wf^- (5 - n) 

In this case, the sum rule becomes: 

If = j (0?l* S (fc)| 2 ^-^), (5-12) 

where } i> B (k) is the Fourier transform of the boundstate wavefunction. The k- integral is only over 
states below the Fermi surface. This makes it clear that 7/ vanishes as v — > 0, and approaches 1 as 
v 1. Roughly speaking, when the density is small, there is a negligible probability of an electron 
being near the origin when the core potential is switched on so the overlap of the unperturbed 
groundstate with any state in which the boundstate is occupied goes to 0. In the opposite limit 
of high density there is probability near 1 of an electron being near the origin. It is interesting 
to consider how rapidly If vanishes as v — + 0. If we assume that the dispersion relation and core 
potential are spherically symmetry, then we may classify the boundstate by its principle angular 
momentum quantum number, L. At small k, ^/ B (k) oc k L , so 

I f oc kf +d oc S+ 2L / d . (5.13) 

In a one- dimensional translationally invariant system with reflection symmetry the boundstate 
can be classified as being an even or odd function of x. If oc v for an even boundstate or u 3 for an 
odd boundstate. When the potential is at the end of a one-dimensional chain with a free boundary 
condition, 

i-i 

* s (fc) = y^smkjip? . (5.14) 

j=i 

Thus * s (fc) oc k as k -> 0, so If OC V . 

The behaviour of the FES exponent, a at v — > follows, in some cases, from Levinson's theorem, 
which determines the behaviour of the phase shift as k — > 0. For an s-wave boundstate, or for a 
one-dimensional problem with the impurity at the end of the chain, the phase shift appro ac hes t t 



at th e bot tom of the band when there is a boundstate. Thus it follows, from Eqs. (2.17), (2.20) 



and (2.21) that aj — » and a e — > 1. Thus If(u>) becomes a step function near its threshold and 
I e (u>) becomes a 5-function. [a = 1 corresponds to a constant Green's function in the time-domain 
whose Fourier transform gives a 5-function.] Thus I e (u), v) approaches, in the v — > limit, the result 
for the empty system: a 5-function of unit intensity. This follows since the groundstate with no 
electrons is the same with or without the core potential; it is simply the vacuum state. Thus we 
expect 

F e (u) 1, (y -> 0). (5.15) 

The step in If(u,v) corresponds to electrons from the continuum falling into the boundstate after 
it is created. The probability for this process, and hence Ff(v) should vanish as v — > 0. 

We have been unable to understand, from general arguments, how Ff(y) approaches or how 
Fe.{y) approaches 1 as v — *• 0. More generally, we would like to know how the functions I e (ui) and 
If(uj) behave even away from the threshold as v — * 0. In particular, it is interesting to ask whether 
there might be some sort of universal scaling form in that limit. This behaviour is only weakly 
constrained by the above sum rules. 

It is interesting to investigate these questions numerically. For a finite system, I(lo, v) is a sum of 
5-functions corresponding to the discrete finite size spectrum. For a large system the intensities of 
the first few peaks, with the boundstate filled or empty, will all be proportional to Ff(v) or F e (y) 
respectively. This follows immediately from the conformal transformation of Sec. II. Thus we may 
conveniently determined Ff i£ (v) numerically, for a large finite system, from the groundstate overlap 
(Anderson orthogonality calculation). That is: 

(0|0)| 2 =Fm(^) 2X . (5.16) 

The resulting functions, Ff >e (v) are plotted in Fig. (|^) for the special tight-binding model considered 
in the previous two sections, in the limit t/V — ► 0. These are obtained from the intercepts of the 
curves in Fig|| As can be seen from Fig.^| at small v, 
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The former equation suggests the scaling hypothesis: 

If{u,v) -> i/ 2 /[(w-w )H, 
for some scaling function /. This is consistent with the sum rule since: 



,oo ,oo 

/ du)I(u>) — v I dxf(x) oc i/ 



(5.17) 



(5.18) 



(5.19) 



In the v — > limit, the threshold frequency c^o approaches e_g — eo, the binding energy measured 
from the bottom of the band. The function f(x) must approach a constant at the threshold, x — » 0, 
consistent with the behaviour of the FES exponent aj — > as v — > 0. At large frequencies the 
function f must vanish sufficiently rapidly for the integral to converge. Apparently the scale over 
which varies is set by v oc vf- If this is also the relevant energy scale for a translationally 

invariant one-dimensional system, or a higher dimensional s-wave boundstate, then we would expect 
the behaviour: 

I(uj,v f ) oc f[(u> - loo)/vf], 

consistent with the sum rule. That is, the threshold peak has a fixed amplitude, but the energy 
scale over which 7(w) decreases scales to as v — > 0. [We can speculate that in our model, due to 
the suppression of the wave function near the end of the chain, the extra factor of v appears in the 
expression for //, which would be absent in the s-wave channel of a higher-dimension system, where 



(5.20) 
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the sum rule should yield v instead of v 3, .] It would be interesting to investigate this behaviour in 
more general models. 

For almost all values of v we found that Ff tB are well described, for our special model, by the 
functions 



,2(1- V* 7 ). 



f» = (1 



\2(1- 



as shown in Figs.^,^. (In two dimensions the dependence rrfi/(V) ~ v 2 ^ l ~ v - 1 could be expected 
from calculations based on the linked cluster approximationEa, but w ith a — 2, not 1/3.) At very 
low densities F e {v) is still well described by (^2l|), as evident from (5.17) and FigM while Fj{y) 
approaches the 20^ 2 -dependence of (5.17) (see FigM) . 



(5.21) 



VI. CONCLUSIONS 



We have investigated the Anderson orthogonality catastrophe and Fermi edge singularity in pho- 
toemission spectrum in a tractable ID tight-binding model of spinless electrons, in the case where 
the core potential produces a boundstate. 

We have confirmed the predicted relation between the scattering phase on the Fermi surface, 
8f, and the Anderson and FES exponents. We have calculated the ratios of intensities of discrete 
adsorption peaks for a finite system, using CFT and checked the formulas at both primary and 
secondary thresholds numerically. We have found that the higher order finite size corrections are 
roughly proportional to 5 F and can be significant in a system as large as several hundred sites. Thus 
they might be observable in a mesoscopic system. The CFT-based relation between the exponents 
and 0(l/Z)-term in the ground state energy shift was confirmed as well. Using the model, we 
obtained the explicit density dependence of the FES amplitude in the whole range of v. 
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Empty 
boundstate 




0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


m = 1 


(S F /n - l) 2 


0.0123 


0.0485 


0.1071 


0.1859 


0.2827 


0.3960 




relative weight 


0.0122 


0.0483 


0.1067 


0.1851 


0.2816 


0.3944 


m = 2 


■(l + (S F /n-lf) 


0.0062 


0.0254 


0.0593 


0.1102 


0.1813 


0.2764 




relative weight 


0.0062 


0.0253 


0.0589 


0.1095 


0.1802 


0.2747 


m = 3 


us f /tt~- + {8 F /n-iyy 

O \ / / V \ - 1 / // 

■(2 + (<5 f /tt-1) 2 ) 


0.0042 


0.0174 


0.0417 


0.0803 


0.1380 


0.2207 




relative weight 


0.0041 


0.0172 


0.0412 


0.0795 


0.1367 


0.2187 


Filled 
boundstate 


V 


0.1 


0.2 


0.3 


0.4 


0.5 


0.6 


m = 1 


(fey*) 2 


0.7906 


0.6079 


0.4525 


0.3236 


0.2193 


0.1374 




relative weight 


0.7868 


0.6051 


0.4505 


0.3222 


0.2184 


0.1369 


m = 2 


i(fe/7r)^ 
■(I + ^/tt) 2 ) 


0.7078 


0.4887 


0..3287 


0.2142 


0.1337 


0.0782 




relative weight 


0.7025 


0.4853 


0.3264 


0.2128 


0.1329 


0.0777 


m = 3 


■(2 + (<5f/tt) 2 ) 


0.6584 


0.4248 


0.2687 


0.1659 


0.0989 


0.0557 




relative weight 


0.6505 


0.4202 


0.2659 


0.1642 


0.0980 


0.0552 



TABLE I. Relative weights of first m excited peaks in the hole propagator compared to CFT predictions of 
Eq.(2.22). The relative weights in the limit I — ► oo are obtained from the best fit to the finite-size values (Fig. 7), 
w(l) = Co + ci/(l — 1) + C2/(l — l) 2 + C3/(l — l) 3 using the method of the least squares to determine Co,.. .,4 from the four 
points with smallest 1/(7 — 1) for each graph. w(oo) = Co. 
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APPENDIX A: FINITE SIZE ENERGY 



In this appendix we wish to demonstrate explicitly the formula relating the finite size correction 
to the groundstate energy difference to the phase shift at the Fermi surface. At the same time we 
will expose a subtlety in the definition of the 0(1/1) term in this energy difference. This groundstate 
energy difference contains a dominant term of O(l). This is non-universal, depending on the ultra- 
violet cut-off in the Dirac fermion theory. It must be subtracted correctly to determine the universal 
0(1/1) correction. To be concrete, we consider a tight-binding chain of I — 1 sites, j — 1, 2, 3, ...I — 1 



When Vj = the single particle eigenstates are sin Ay. A simple way of determining the allowed 
wave-vectors, k, is to imagine adding two "phantom sites" at j = and j — I and then imposing 
the boundary condition: 

= "4>l = 0. 

The Fourier expansion of ipj in terms of creation and annihilation operators then involves sin(fcj) 
with: 

k n = Tvn/l, m = 1, 2, 3, . . . I — 1. 
The groundstate energy for V = is: 



with e(fc) = — 2icosfc. In fact, this discussion doesn't depend on the form of e(k) and can be applied 
immediately to a more general Hamiltonian with longer range hopping provided that the boundary 



condition of Eq. (A.2) applies. The only property of e(fc) that we will use is that its derivative 
vanishes at k — 0. 

In the I — > oo limit, the Fermi wave- vector is 

k F = n lim N/l. 

I — >oo 

For a finite system, there is an ambiguity of 0(1/1) in the definition of /cf since it may be chosen 
anywhere between irN/l and n(N + l)/l. It turns out to be convenient to choose it to lie exactly 

halfway between the N th and I st level: 

k F = ty{N+1/2)/1. 

This gives the model an approximate particle-hole symmetry, in the vicinity of the Fermi surface. 
[Only at half-filling does this particle-hole symmetry become exact.] We regard k F as being held 
fixed as I is varied, for purposes of determining the term of 0(1/1) in the groundstate energy. Thus 
the quantity (N + 1/2) jl must be held fixed. In practice, for numerical simulations, this is not 
particularly more nor less difficult than holding fixed the actual density, N/(l — 1). For instance, to 
obtain k F — tt/4 we may choose the number of sites I — 1 = 4iV + 1 for arbitrary positive integer N. 
The continuum limit Dirac theory is defined by only keeping wave- vectors near ±k F , writing: 



where iph,R are left and right moving Dirac fields. The boundary conditions of Eq. (A2) imply 



^t(0) + ^«(0) = 
- ik " l ^ L {l)+e ikFl i> R {l) = Q. 



Using Eq. (A6) the last equation gives: 

^ L (l) - ip R (l) = 0, 
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corresponding to t he "same" boundary conditions at both ends as discussed in Section II. k F was 
chosen in Eq. (A6) above in order to obtain this boundary condition on the Dirac fermion. 

Including the scattering potential, Vj, the single-particle eigenstates still become asymptotically 
plane waves, sin[fcj + 5(k)] (at distances large compared to the range of V), where S(k) is the phase 
shift. Thus the allowed wave-vectors are now: 

kn = kn 3{kn)/l] 



with k n = Tin/ 1. It is important to note that the argument of S in Eq. ( A10 ) is k n , not k n . Thus, 
to 0(l/l 2 ), we may write: 

k n = f{k n ) =k n - 8(k n )/l + 5'(k n )5(k n )/l 2 . 

The groundstate energy can thus be written: 



E = 



This can be evaluated using the Euler-MacLaurin expansion: 



J] F(rn - 1/2) = f 
i Jo 

m— 1 u 



dxF(x)--[F'(N)-F'(0)]+O(F"). 



Setting: 



F{n- 1/2) =e[/(7m/0], 



where the function, / is given by Eq. (All), we obtain the convenient result, 

F(N) = e[f(k F )]. 

Thus, to 0(1/0: 



E = f 
Jo 



where 



dne{f[Ti(n + 1/2)//]} - v F ir/(24I), 



v F = e'(kF), 



and corrections of 0(l/l 2 ) have been dropped. Now it is convenient to change integration variables 
to: 



giving: 



k = 7r(n + l/2)/l, 



Here the lower limit of integration has been shifted by n/2l. This is justi fied since e(k) is quadratic 
at k — > 0, producing only corrections of 0(l/l 2 ) to Eq. Using Eq. (All) and expanding to 0(1/1) 
we obtain: 



E = l 



dk 



<k) _ e'(k)5(k) + e"(k)S 2 (k) + e'(k)6'(k)5(k) 



I ' 2l 2 

Integrating by parts, and using e'(0) = 0, we finally obtain 



~2M 



dk 



E = I I —e(k) - - I de5(e) + ^ {- 



TYVf I 1 



S(k F ) 



24 



l 2 



Here eo = e(0) and £f = e(fe_p). The first term, of 0(1) is the bulk groundstate energy. The 
second term of O(l) is a well-known result referred to as Fumi's theorem. Note that these terms 
depend on e and 5 over the whole band. On the other hand, the final term, of 0(1/1) depends 
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only on data right at the Fermi surface, namely vf and S(kp). The 1/24 term is the well-known 
conformal field theory result for open boundary conditions. The additional {S(kp)/iv] 2 term gives 
the effect on the groundstate energy of changing the boundary conditions. This formula can be 
checked for the nearest neighbour model, with e(k) = 2tcos(fc), and free boundary conditions. The 
exact groundstate energy is given by a geometric series: 



E = t ■ 



t sin kF 
sm(%/2l) 



I irv F 

— vf +t — — 

7T 24/ 



o{\/r 



with vp = 2t sin kp ■ 

We note that adopting a different definition of kp, such as tvN/(1 — 1) in the one dimensional 
model, corresponds to ad ding a fc- independent term to the phase shift [(u — 1/2), where v is the 
density, in this case]. Eq. ( A22) still holds, when written in terms of this redefined phase shift. This 
discussion was given for the case of a one dimensional tight-binding model with free ends but it can 
be easily generalized, for example to the s-wave sector of a three dimensional spherically symmetric 
continuum model with a vanishing boundary condition on the surface of a sphere if radius, I. With 
an ap propriate definition of kp [essentially defining, to 0(1/ 1) what is held fixed as I is varied], Eq. 
(A22) is again obtained. 



We can easily generalize Eq. (A22) to calculate the energy of a state with n extra electrons added, 
with n held fixed as N and I — > oo. This gives: 

n 

E n =E + Y^ #f - 5{k F )/l + 7r(m - 1/2)//] 

m— 1 

71 

= E + ne F + {vpTv/l) ^ (m - 1/2 - 5{k F )M + 0(1/Z 2 ). 



(A22) 



(A23) 



Hence 



E n = / 



dk . 1 [ ep , . -kvf 1 

— e(k) / ded{e) + ne F + —— { - 

TT TV I Z 



5(k f 
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+ o 



(?) 



(A24) 



APPENDIX B: DISPERSION LAW AND WAVE FUNCTIONS IN A ID TIGHT-BINDING CHAIN 



Here we calculate exactly the phase shift, finite size spectrum, eigenstates and overlap s fo r the 
one-dimensional tight-binding model with an impurity potential at one end, given in Eq. ( |3.l[ ). 
It can be easily seen that the eigenstates can be written exactly in the form: 

tyj oc sinfc(j — I). 

This wavefunction trivially satisfies the lattice Schroedinger equation (for arbitrary k) at all sites 
2,3, . . ., l-l, with: 

e(k) = —2t cos k. 

The Schroedinger equation for the first site determines the allowed values of k: 

-£* 2 - V*i = 6*1. 



Inserting Eq. (Bl) and (B2l), we obtain: 



sin k(l — 1) t 
sin kl V 

For sufficiently large V/t there is also a boundstate with wavefunction: 

Xj oc sinhft(/ - j), 

where k > in order that the wavefunction decreases with increasing j. This has energy: 



(Bl) 

(B2) 
(B3) 

(B4) 
(B5) 
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e_B = — 2i cosh re. 

Again the Schroedinger equation is satisfied automatically (for any re) at all sites except the first 
which gives the condition determining re: 

sinh k(1 — 1) t 
sinh nl V 

For Z >> 1, this gives: 

e _K = t/V. 

Since re must be positive, we see that there is only a boundstate solution for V > t. In this case: 

e B = -{V + t 2 /V). 

For the continuum states we may calculate the exact phase shift, defined by the form of the wave- 
function for I » 1: 

*j oc sm[kj + 6(k)]. 



From Eq. (Bl) we see that: 



-kl = S(k) + nn, 



for some integer n, in the limit I — > oo. Substituting this into Eq. (B4), we obtain: 

sin[<J(fc) +k] _ t 



sin[<J(fe)] V 



This gives: 



8 = arctan 



sin k 



t/V — cos k 



We now consider in more detail the spectrum of the finite system, of I— 1 sites. With no potential, 
there are I — 1 band wave-functions with: 

k n = rm/l, n = 1, 2, 3, . . . (I — 1). 



Including the attractive potential, we see from Eq. ( B7) that there is a solution of the form sinh re(j — 
I) for t/V < 1 - 1/7. We label this solution h and f 1 . For this range of t/V Eq. (§) has 
only 1 — 2 solutions, fo, &3, . . . In particular, for V/t — > oo, these / — 2 solutions become: 

k n = (n — l)vr/(/ — 1), corresponding to a chain with free boundary condition at both ends and I — 2 
sites. 

The normalization of the band states can be calculated exactly in terms of k using: 



^sin 2 fc(j - I) 



3=1 



sin k(l — 1) cos fcZ 



sin k 



Similarly the normalization of the boundstate is determined by: 



^2 sinh2 K 0' - 0=o 



3=1 



sinh re(Z — 1) cosh kI 



sinh re 



The overlaps of band wavefunctions with and without the potential can be calculated similarly 
using: 



S ' sin fc(j — I) sin fc(j — I) 



3=1 



sin fc(i — 1) sin Zfc — sin Ik sin(Z — l)fc 
2 [cos — cos k] 



If and k are allowed wave-vectors corresponding to potentials Vi and V2 respectively, then using 
Eq. (pi), 
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y sin k(j I) sin k(j - I) = (^-Vi)sinfc(f-l)sinfc(i-l) (Blg) 
^— ^ 2t[cos fc — cos fc] 



Thus we obtain the extremely useful result: 



<*\*>=y^m£M, ( Bi9) 

e(fc) - e(fc) 

where e(fc) = — 2tcosfc is the band energy and 

c(k)= - ^infc(z-i) (B2Q) 

•^/(Z — 1) — sin fc(i — 1) cos klj sin fc 

The corresponding result involving the boundstate follow immediately upon replacing k by in. e(k) 
simply gets replaced by es = — 2icoshK and C(k) by 

C B ee V2sinh K (Z-l) (B21) 

v/ sinh k(1 — 1) cosh kI / sinh k — (i — 1) 

To calculate the overlap of the V = groundstate with an arbitrary state with V 7^ 0, we simply 
set Vi = and V2 ~ V in the above formula. This remarkably simply form for the overlaps of single 
particle wavefunctions leads to enormous simplification in the calculation of the overlap of the Bloch 
determinant multi-particle states, as shown in Section IV. 
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